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ABSTRACT 


Woaverrebacea velocity ftlue@tuattens over ocean waves and 
wave heights as measured from FLIP during BOMEX are examined 
using phase-amplitude results which are based on joint prob- 
ability density function-conditional mean function (JPDF- 
CMF) analyses. Results are compared with predictions from 
various wind-wave coupling models. 

Legulesmare examined in detail and consistent departures 
meamechneOry are noted. An attempt is made to qualitatively 
M@ewermine the effect on specific results of the moving, but 
mematively Stable, sensor platform, FLIP. 

It 1s concluded that the interaction between the wave- 
menecd Motion and airflow turbulence had a significant 
effect on the observed wave-related fluctuations. The 
SeEccts of FLIP on the results aApPeakee oO De Milita lyon 


@hese results. 
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I. INTRODUCTION 


The atmospheric boundary layer above land has been 
aiojrougiiy studied and 1s well known, However, over the 
oceans, where the mobility of the underlying surface has to 
be considered, the investigation and knowledge of the 
mermueture and dynamics of this layer is still incomplete. 
Sontlicts have existed in the theories describing the 
basic dynamic interaction between surface waves and the 
airflow above them. 

iheoretical models, with’ valid, comparable observational 
meewlts, are required to provide detailed descriptions of 
the structure and energy transfer mechanisms of the near 
Stirface layer over the ocean. Examples of geophysical 
problems in need of this information include: long-range 
memospherie prediction models which require, for the nea 
aimrace layer, a better specified input of the kinetic 
energy losses and the thermal energy gains and surface 
wave predictions requiring wind speeds at specified levels 
irethne surface layer. 

The primary purpose of this study is to perform joint 
probability analyses on velocity data obtained in the near 
Surface layer and on simultaneous wave data. vom these 
analyses, interpretations will be made on the phase- 
amplitudesrelationships, nonlinear properties, and some 


determination of the influence on turbulence present. These 








jeoults are compared to existing wind-wave coupling theories 
and to other observational results which have only recently 
compared favorably with these theories. 

PPC nent: proplem im tie collection=of meteorological 
@eea in the near surface layer over water waves in the ocean 
mmeene etfect of the motions of the sensor platform. The 
data used in this investigation were accumulated using 
S-7wsoors mounted on a platform which, although relatively 
Stable, did experience some motion in response to the waves. 

A parallel purpose of this study will be to attempt to 
@eamermiine dualitatively, if not quantitatively, the effects 


feeplattoOrm motion on the results from these specific analyses. 





II. BACKGROUND 


A. SEA-AIR INTERACTION EXPERIMENT 

In May of 1969, the "Barbados Oceanographic and 
Meteorological Experiment" (BOMEX) took place in the 
Pralemipiic Ocean egeae | Or Ene ast of Barbados, West Indes. 
BOMEX was a joint investigation by several agencies and 
universities, including a team from the University of 
fienigan which collected the data used in this study. The 
overall purpose of BOMEX was to establish the energy budget 
of a volume 500 km square and 5 km above the ocean surface 
to 500 m below it. 

As one of the largest environmental research programs 
ever conducted, BOMEX contained over 80 subprograms. Of 
articular interest to this study was the “sea-air interac- 
tion'"' subprogram. One aspect of this subprogram was to 
collect data in the near eri ace layer. Various oceano- 
graphic research platforms were used in these measurements. 
The University of Michigan measurements were taken from the 
Scripps Beene on of Oceanography Floating Instrument 
Platform (FLIP) which was located near the center of the 
volume at approximately 16°N latitude and 57°W longitude. 
Bronson and Glasten (1968) have described FLIP in some 


detail. 


B. THEORETICAL BACKGROUND 


Theoretical investigations of wind-wave coupling have 


progressed from the quasi-laminar model proposed by Miles 





aS 





Mee) to fhe recent tumbulence models, such as formulated 
aeceimoy (1970). ithe close agreement between recent 
Seeemvatitonal results and numerical solutions from recent 
wind-wave coupling models by Yefimov (1970), Reynolds (1968), 
amd several others provided the primary motivation for the 
@eeaen of this study. 

Most pemeoreticalwsfoumulations of linea, theory for wind- 
wave coupling have several common assumptions such as two- 
dimensional flow and a specified mean wind dependent only on 
Meieht, and other common simplifications such asa single- 
Somponent wave field and suitable linearized forms for the 
aomiming Equations. A distinguishing feature of the 
formulations by Yefimov and Reynolds, and more recently 
Myerotautfer (1973), has been the allowance for the interac- 
meonmepe tween the wave-induced motion and the turbulence in 
mmemoverlying Shear flow. In contrast to this, Miles (1957) 
Mermected turbulence except for its role in determining a 
assumed logarithmic mean profile. Davis (1972) provides 
description and evaluation of the approaches in the more 
imeeent turbulence mocionee 

Hitless’ theory 1S a quasi-laminar wave generating 
mechanism wherein the wave-induced motion is assumed to 
mollow the potential flow theory. The existence of the 
Sebeiealelevel, which corresponds to the level where the 
wind speed is equal to the phase speed of the perturbing 
wove tomelemkey tOuMales theory, Below the critical 


level the airflow would appear to be opposite that of the 


Jen 





mares propagation Since the wind speed is less than the phase 
speed of the wave. Above the critical level the wind speed 
Beveeceds the wave speed; hence the flow 1S opposite to that 
Below DUt in the same direction as the wave. Miles’ theory 
is based on the premise that waves only modulate the mean 
wind flow over them and that turbulence is not influenced; 
mitsethe Reynolds stresses are only dependent on the height 
apeme the mean wave surface. Figure 1 illustrates the 

Pmase and amplitude features of two components of the velocity 
and the stress uw, with respect On tle waver as prealeted 
Byepetential flow theory. These predictions apply to 
meueGtuations observed at a fixed level in a shear flow 

above a progressive wave. 

In Yyetimov's formulation the allowance for turbulence 
in the airflow rr ae accounting of the nonlinear 
Berms 1n pepcquecions Cf motion=. [ne model provides stnat 
mmoamics at the critical level are not the only reason for 
hemes related motion to deviate from potential flow predic- 
tions. Several Gop etiees Mie vetimNoy saSOlutLlons are due 
solely LOuEhCuMecscribeCd Interaction between the prescribed 
motion and SBM Oe: AS Suekhye they occur independently 
Seewemne height or of the existence of the critical level. 
nese features were the variations with height of the phase 


differences between the velocity components and the wave. 


C. OBSERVATIONAL BACKGROUND 
Theoretical investigations of wind-wave coupling have 


seemingly progressed at a steady rate; however, until 
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Figure 1. Phase-Amplitude Relations Predicted by Potential 
Flow Theory for the Two Components of Velocity 
(u and w) and the Stress uw, with Respect to the 
havem@inonm Wavtason and Prank, 1973). 
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recently the availability of overwater measurements has 
Decmeimsurticient tO C€ither verify or modify existing 
theories. Over the last six years the situation has 
improved considerably. Observational studies of measure- 
ments over natural waves in the near surface layer have 
Deen described by Davidson (1970), Pond, et al. (1971), 
Memeo et al. (1972), Davidson and Frank (1973), Thompson 
oer 2) and a number of others. Of particular interest to 
this thesis are the investigations of Davidson and Frank, 
and Thompson. Their studies are similar in types of analysis 
midmreatures identified to those used in this study. 
Davidson (1970) ae Simultaneous measurements of 
velocity fluctuations and wave heights at two separate 
Mewes On Lake Michigan. These data were analyzed by 
fweetral methods and showed rather clearly that the air- 
flow was being influenced by the underlying wave field. 
Davidson and Frank (1973) re-examined the Lake Michigan da a 
using joint probability density-conditional mean aualysis 
mmgocedures. the earlier spectral result Swere “compared art 
the wind-wave coupling predictions of Miles' theory while 
the JPDF-CMF results considered the role of turbulence in 
these same data. Recently the results of these analyses 
were used for comparison in a numerical turbulence model 
by Stauffer (1973). The model used by Stauffer is very 
Simla to Yefimov's. The comparison of the velocity 
components, u and w, and the stress uw, showed agreement 


between the predicted values and those observed. 


14 





Thompson (1972) used essentially the same analysis 
procedures as those used by Davidson and Frank. In his 
analysis he used some periods of the same data analyzed 
Momenis investigation. Thompson interpreted JPDF-CMF 
based phase-amplitude results with respect to potential 
Pew predictions and the role of turbulence in wind-wave 
eempling, He observed that the wave-related motions in the 
Mear surface layer were essentially nonpotential. 

iigese two investigations are of primary interest in 
memwect to this investigation for several reasons. Both 
iimestigations used analysis procedures which were the basis 
memetnose used in this study and therefore result in the 
fale Parameters for comparison. The Lake Michigan results, 
Memrnaich these results can be compared, are in agreement 
Peem those predicted by a numerical model. In addition, 
the Lake Michigan data were accumulated from a fixed 
mieenument platform. Comparison of the results from this 
investigation with those from Lake Michigan should provide 
some indications as to the effects of a moving instrument 
platform. And finally, Thompson, in his analysis of BOMEX 
data, made minimal consideration for the DOssiole Iieluence 
GreELIP'ts motion on the results. 

Thompson followed the analysis by Rudnick (1967) in 
PeoteCrligethie possible influence of platform motion on the 
results. Thompson indicated, on the basis of Rudnick's study, 
that his observed fluctuations in u and w were approximately 


a” times larger than could be induced by platform motion. 


PS 





fevever, Pond (1968) motes that there are marked differences 
between measurements from buoys and fixed masts. Pond 
states that provided some damping (by design or mooring) is 
done of the vertical motion and provided that the amplitude 
On the Sseriiarmestrit is smaller than 10°-20°, the effects 
Gmpnoy motion are fairly smail. He shows that if these 
Semadtions are met the effects of buoy motion are "probably" . 
10-20%. 

Provided? Pie ppemaixww 1S, ad fadvrly simple attempt 
pemence representation of the apparent wind field in rela- 
tion to a moving spar buoy. An attempt will be made in 
the presentation of results to use these representations, 
and the results of Rudnick and Pond to determine qualitative- 


ieeene effect of FLIP's motion. 


Ihe descriptions in Appendix A were provided through personal 


Scr Uitavronmomancameorrespondence with Dr. A. W. Green of the 


Department of Meteorology and Oceanography, University of 
Michigan. 
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Lite Usa COLU EC ON 


University of Michigan personnel obtained 54 hours of 
data in 40 separate observational periods during the BOMEX 
peopemiment, Detailed descriptions of the sensors, sensor 
arrangement and mounting, and recording equipment are 
memlaple im reports and studies by Portman, et al. (19/70), 
aaeinompson (1972). This information will only be briefly 
mentioned in the following paragraphs. An aspect of signifi- 
Sime importance for this study is consideration of FLIP and 


Mmeceise as an instrument platform. 


A. INSTRUMENTATION 

ete Lullteeconstant temperature anemometer systems 
Pemenmounted at the 2, 3, 6, and 8 meter levels on the 
Memueical mast (See Figure 2). These sensors were capable 
of measuring the three fluctuating velocity components 
(u, v, and w) at each level. 

A resistance wave gage, provided by Dr. R. E. Davis, 
Scripps Institution of Oceanography, was used to measure 
wave heights. The wave gage was positioned about five 


meters inboard of the vertical instrument mast. 


pe FLOATING INSTRUMENT PLATFORM (FLIP) 
In the previous section the motions of the spar buoy, 
wichedsetlllP were discussed from a general, theoretical 


standpoint. The following paragraphs are a more specific 
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Bure 2. View of FLIP Showing the Arrangement of 
Somaors during BOMEX. 
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Figure 3. Top view of FLIP Showing the Mooring 
Procedure Employed during BOMEX. 
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fesepiption of BLIP and its use as an instrument platform 
during BOMEX. 

Mollo-Christensen (1968) reported that a minimum of 
fbh> tortion to the natural wind flow by FLIP's super-structure 
Peould occur if the deck of FLIP were positioned into the 
a. This was determined in wind tunnel tests using a 
Secale model of FLIP. Mollo-Christensen recommended that 
mites vertical Sénsor mast be positioned about 15 meters 
Saepoard of the hull. 

Figure 3 shows the method utilized during BOMEX to 
em in FLIP's position with respect to the wind direction. 
Mime tug waS positioned about 800 meters downwind of FLIP. 
The tug was unable to maintain constant tension on the 
memineg cable and this has been discussed by Superior (1969). 
meead result of a varying tension on the towing cable and the 
meotrcioning method in general, axial and lateral motions of 
Wiper occurred. 

Untonrtunately, real time measurements of FPLIP's motion, 
due to the features described above and those caused by 
wave mo con’, Bee NOtwaveiiabre., because Of this Tt 1s 
doubtful that the exact effects of FLIP's motion on the data 
seulected can be accurately described. 

Another problem of concern to data interpretation was 
that provisions were not available during BOMEX for adjusting 
the vertical position of the sensor mast, independent of 
Orton well seecaused the sensors to be out of level 


dite potrelOmmor seach peTlod . 
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IV. THE PREPARATION AND ANALYSIS OF DATA 


Joint probability density--conditional mean (JPDF-CMF) 
agmalyses were uSed to identify features in the fluctuating 
velocities which were associated with the waves' presence. 
Phase-amplitude results were then based upon the JPDF-CMF 
analysis. 

In general, the procedures used in the analysis of the 
data for this study were the same as those originally described 
Bbyetiolland (1973) and applied to these kind of data by 
Thompson (1972). Improvements were made in the procedures 
for defining phase-amplitude information from the JPDF-CMF 
aomeeiyses. ihe resolution of the phase-amplitude results 
was increased to enhance the shape and phase relationships 


mieene curves. 


A. INITIAL DATA PREPARATION 

Several basic data processing steps! preceded the 
awoeredt10On of the JPDF-CMF analysis procedures. These 
maemnded digitizing, application of an analog filter, time 
@emrections, scaling of data to engineering units, and 
meveral other procedures described in adequate detail by 
pumghiam (1972). An iteration of the details of these 


peeliminary steps will not be included in this report. 


initial data preparation procedures were developed by the 
Depanunene Of Meteorology and Oceanography, University of 
Mictiecai mech the —dtrection of Professor D, J. Portman. 
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B. BAND PASS FILTERING AND NORMALIZATION 

In order to remove high and low frequency contamination 
mee to Locus the analysis on fluctiations near the frequencv 
band corresponding to the water waves, a band pass numerical 
meeerse transtorm type filter was applied to all data. This 
ype Of filter could be very useful in also removing low 
imeauency Contamination to the data from FLIP's motion. 
The filter was designed on the basis of available spectral 
meomles and had a frequency band with 0.1 and 0.8 Hz as 
the low and high cut-off frequencies, respectively. The 
imerter’'s response curve appears in Figure 4, 

Prom tO application of the JPDF-CMF procedures the 
data were normalized by dividing the deviation from the 


mean by the standard deviation. 


10 


J n (Hz) ! 10 


Figure 4. Response Curve of the Band Pass Numerical 
hivercewiransrorm filter. 
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tee Pir-CMF PROCEDURES 

fMespamtlcilar procedures used in the JPDP-CMP computation 
were developed from descriptions by Holland (1973). A brief 
beemcomplete discussion of these procedures 1s given by 
Beavadson and Frank (1973) and Thompson (1972). The JPDF-CMF 
eeeiational DEOCcdutecmeaiicmeanplied, CG three variables , 
meee im the joint probability density function (JPDF) for 
ame, Of variables with the conditional mean of a third 
wale as a function (CMF) of the first two variables. 
Representation of the resulting trivariate statistical 
Peewerouship and further discussion of the JPDF-CMF analysis 
meemeontained ian later sections. The JPDF-CMF computations 
were made from the normalized values of the original varia- 
Muecwoampled at five points per second, with a majority of 
micemmecords consisting of 6062 points or approximately 20 


mmmees of data, 


Weeeeitaok-AMPLITUDE COMPUTATIONAL PROCEDURES 

Hierle various features of the fluctuations in the airflow 
above water waves can be identified by JPDF-CMF analysis, 
meoinpeitication of the method of display is desirable in 
Greer tO concentrate attention on the essential variations 
and eliminate redundant background information. For this 
Saray, the emphasis is on the statistical dependence of 
turbulent variables on the phase of some reference variable. 
Based on JPDF-CMF procedures, Holland (1973) described a 


method to determine this phase-amplitude information. 
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Pollar coordinates representing the amplitude and phase 
angle of a variable such as the wave height (n} can be 
determined from the JPDF array for the variable, n, and its 
time derivative, nj. In amplification, any conditional mean 
function of n and 4 is a CMF of the phase and amplitude 
of the n fluctuations because for each value of the dependent 
variable represented by the CMF, the coordinates in n, n 
Space are determined by the amplitude of the n fluctuation 
Peering at that time and by the phase angle within that 
Mmemeedation, As an example, the conditional mean function 
w(n,n) would illustrate the dependence of the» vertical com- 
ponent of wind velocity on the phase and amplitude of the 
wave. | 

An example of the JPDF-CMF array and the resulting phase- 
mipgeeuide Curves appear in Figures 5 and 6. Indicated on 
the array are the numbered 30° segments and the class limit> 
Wmrermwere used to approximate the 1.50 circle. The small 
popiede class (amplitude < 1.50) contains about 3/4 of tI: 
observations and the Goce walpleneice «class —(amplicude > vies) 
@emeains about 1/4. Selection of these class limits is some- 
Miteeearbitrary. The phase is defined in terms of the 30° 
segments measured counterclockwise in the Fron Gooremace 
Syetem (JPDF-CMF array) from the positive 2 aoaiice 

Using Figure 5 as an example, the following procedures 
Hewes Used £o Obtain the resulting phase-amplitude curves. 

Peipbewcreroures = (IPD array) in each joint class inter- 


val are summed over each 30° segment and amplitude class. The 
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iemmcedensitiés in joint class intervals along the segment 
ingiessare divided by percentages and are added into the 
Seanemts Which share these joint intervals. Each probability 
Poimers then divided by the sum over the entire amplitude 
Glass. 

me Products of the dependent variable (CMF) and 
DPme@ptbality density (JPDF) in a]] joint intervals of each 
30° segment and amplitude class are summed and the sum 1s 
divided by the total probability in that segment and 
amplitude class. 

Fesults Of these procedures yield a single curve for 
each dependent variable, showing it as a conditional mean 
Mmiterron Of the phase angle of the reference variable, which 
feoteoo illustrated as a single curve (Figure 6). 

leeviously, investigators using these procedures 
(Davidson and Frank, 1973, and iionpson. 2972) pertormed “the 
Summations over octants (45° segments). In this study, using 
30° segments, a considerable improvement over the previous 
studies in the resolution of the phase-amplitude curves was 
achieved. Phase-amplitude results using octants are shown 
waericure 7 for comparison with those in Figure 6, 

Additional aspects of the phase-amplitude results are 


presented during a discussion of results. 
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Pee 6. Phase-Amplitude Results Showing the Resolution 
Achieved by Summing over 30° Segments. 
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Figure 7. Phase-Amplitude Results Showing the Resolution 
Achieved by Summing over 45° Segments as used 
by Previous Investigators. 
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Vo PResrNt VON GGEREOUIOES 


The data from 14 periods were analyzed during this 
mimectigation. Ihey were collected on five separate days 
and represent over 15 hours of observations during BOMEX. 
The periods were from 34 minutes to 80 minutes in length. 
mmabyses were performed on data from the three, six and 
eight-meter levels with one period consisting of simultane- 
@ecedata from both three and eight meters. For the purposes 
of JPDF-CMF and phase-amplitude analyses, each period was 
Seevided into subsets (2 to 4), each of which was 20 minutes 
@aeless an Length. The last subset in each period was 
mMeulially less than 20 minutes long. Table I gives a summary 
merall periods analyzed, indicating times of observation, 
Mement analyzed, mean wind speed at that height, and several 
Sener parameters. 

General meteorological conditions, wave heights, and 
@ueresponding times as reported by University of Michigan 
personnel during BOMEX are given in Figure 8. The general 
Memes pneric conditions during all periods could be described 
as stable or near neutral. Profiles of wind and temperature 
Measurements from two meters to 16 meters can be found in a 
Study by Superior (1969). There were no periods of active 
Wave growth or rapid decay. The conditions for the area in 
which FLIP was located could be described as a region of 
light to moderate winds into which swell was propagating. 
Thus it was fairly representative of the surface layer over 


Bhemoeeonmana relevant, tO this discussion. 
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eieeresults ot the analyses of the period numbers marked 
tameiean asterisk in Table I will be those interpreted in 
@etart in the following paragraphs, and include period 
fmmeers 1, 6, l2Za, 12b, and 14. Results from these particu- 
lar periods were chosen for several reasons. Periods 1 and 6 
Meemesent very nearly the extremes in observed mean wind 
Speeds, while periods 6 and 12 represent the extremes in 
Seserved wave heights. Periods 12a and 12b, 8 and 3 meter 
Meamete results, respectively, were also chosen because they 
were the only multiple level measurements available and, 
Memec, can bevinterpreted with respect to phase-amplitude 
@mmmecs with height. They also should illustrate to a 
Mmeemmer extent the effects of FLIP. Period 14 data was 
oommeeceted at the 3 meter level and had winds and wave 
iememes Detween the extremes. This period corresponds to 
Beraod 3 ener by Thompson (1972) in which the motions 
of FLIP were essentially neglected but herein will be re- 
exomumed with consideration of sensor platform motion. 

Data sou sidlened In tiomammn~ses Were. the U, Vv; and w 
velocity components, the stress uw, and the wave heights. 
Associated en each of these variables is a parameter To 
mee 15 defined as the "statistical eddy period'' and is 


@eeermined by the ratio 


Bata Zan 7.2 


where o and o are the standard deviations of the variable and 


Pmecmivearive=. Am indicator of the appropriateness of the 
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ailysisS 1S to compare Ten POReuNeracpenuenl Variable 
Wee ectcrye with that ot the™reference variable, the wave 
ements If the Statistical eddy periods differ signifi- 
eemitly, the phase-amplitude information for the reference 
meaiable has a greater possibility of only representing the 
een eo: numerous random fluctuations. Statistical eddy 
meeods are listed on each of the figures depicting the results 
pomeeros in fact, attest to the validity of the analysis. fhe 
Seataistical eddy period for the product uw is less than that 
Memetie wave because uw in a linear system, at least, would 
appear as a second harmonic. For this reason, we 1S Oe 
ested for uw. 

meceptance 15 Piven, in the interpretations of the phase- 
amplitude results, to Kendall's (1970) Suggestion that non- 
sinusoidal wave forms are evidence of non-linearities anes, 
Meee a non-negligible role for turbulence. It should be 
Mmoeea that the "critical" level was found to be above 10 
meters for all periods. In the following paragraphs descrip- 


tions of the periods discussed above are given in detail. 


eee RIOD 12 

hits period 2s Unaque in that it is the only period 
available from the data examined that contained simultaneous 
information from two levels and, as such, accounts for a 
major portion of this section. 

The duration of the period was about 80 minutes and, 


hence, four subsets are available for each level. Wave 


op 





fements and wind speeds were observed (Figure 8) to have 
peeneincreasing in the four to five hours preceding this 
period and reached a peak just prior to or in the early 
tometen of the period. ~The wave heights for the period 
tease observed to have been approximately seven to 11 feet 
Wienemean wind speed of 10.0 M/sec at the 8 meter level and 
9.8 M/sec at 3 meters. 

Figures 9 through 12 contain the phase-amplitude results 
Mamie S meter level while Figures 13 through 16 depict 
the 3 meter level. Only the phase-amplitude relationships 
Gumme W, and the wave will be discussed in detail. 

the u components at both levels and for all subsets 
femerae the period are surprisingly similar. The phase of 
Beemcomponent relative to the wave in all cases is as 
memmmcam pe predicted by potential flow theory except that the 
oeeeMoswot uw have shifted back slightly. Potential flow 
MredtetsS, as shown in Figure 1, that u and the wave shoul 
be exactly 180° out of phase with each other. The phase 
Pim@eeobserved in’ these results shows u to be trailing th 
Sucsemor trough of the wave by about 15° to 30°. This shift 
was observed by both Davidson and Frank (1973) and Thompson 
(1972) and they associated the shift with momentum transfer. 
itesamplitudes of the u component fluctuations decreased 
slightly from 3 to 8 meters. 

The phase relationships in w are the most interesting 
of the results from this investigation. Between the 3 and 8 


eter wlevelmenere 1s a phase Change for w of between 90° 


oY 





fiwmemres 11 and 15) and 150° (Figures 9 and 13). As can be 
seen for the 8 meter level, the phase relationships of w 
with the wave and u are not consistent with potential flow 
mreory With the possible exception of subset 3 in Figure ll. 
Samerevidence is found to indicate the role of turbulence 
in that the wave-forms of w appear to be non-sinusoidal in 
maiesets 1. 2, and 4 (Figures 9, 10, and 12). The fact that 
the u and w components at the 8 meter level are nearly in 
phase was also noted by Thompson (1972) and was attributed 
to turbulent interaction with the wave-induced motion. In 
mcet 4 (Figure 12) it can be seen that there are no phase 
differences between u and w at the negative extrema; however, 
at the positive extrema w is leading the u by about 30°. 
Wemential flow theory predicts that u and w would be in 
Selacrature aS indicated in Figure 1. 

Non-sinusoidal wave-forms are also evident for the 
Meietwer level as seen in Figures 13, 15, and 16. At this 
Mevel the departure of the vertical fluctuations from 
potential flow are most noticeable. Instead of being in 
quadrature as eee a Umeneaiw are almost 180° out of 
phase. The maximum w fluctuation occurs when u is a minimum. 
This shift in w was also noted by Davidson and Frank (1973) 
mame results of their period 2 at the 1.5 meter level. 
iiey attributed this shift in the w component to “critical 
level" dynamics as a prediction in the quasi-laminar theory. 


However, with a change in the phase of w with height it can 





femassumed that it was caused by turbulent interaction in 
Bimes Case. 

This period provides some evidence that perhaps FLIP's 
Hemelom did not appreciably affect these results. First, 
since FLIP is a muicidebody and the venmtrcal sensor mast as 
ome attached, errors introduced by FLIP at one level 
would also be introduced at the other level. In other 
Mendes. the £luctuations introduced by FLIP at one level 
would also be evident at the other and thus no phase shift 
could occur unless it was present in the airflow. From the 
discussion in Appendix A, FLIP could be considered an ampli- 
meeeiinder certain conditions and, as such, could introduce 
Mime Shifts into the results. However, since the results 
fimemiiccers aré in agreement with the results at a lower 
level from a fixed platform poanaecon aid Clamita camtke a, 11997 5.) ee tetas 
pemomertlity will not be considered. FLIP may introduce some 
SGmpor into the magnitudes of the various fluctuating variables. 
However, it is doubted that that error could be determined 


mEomeeernese results. 
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Figure 9. Phase-Amplitude Results for Period 12a, Subset 1 
(S-meter Wevel). 
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Figure 10. Phase-Amplitude Results for Period 12a, Subset 2 
(8 meter level) 
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SMALL LARGE TOTAL 


Figure 11. Phase-Amplitude Results for Period 12a, Subset 3 
(8 meter level). 
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Figure 12. Phase-Amplitude Results for Period 12a, Subset 4 
Pomme tcr elec!) . 








SMALL LARGE TOTAL 


Figure 13. Phase-Amplitude Results for Period 12b, Subset 1 
(3 meter level). | 
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~ SMALL LARGE TOTAL 


Figure 14. Phase-Amplitude Results for Period 12b, Subset 2 
(3 meter level). 
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SMALL LARGE TOTAL 


Pane 15. Phase-Amplitude Results for Period 12b, Subset 3 
(3 meter level). | 
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Pigure 16. Phase-Amplitude Results for Period 12b, Subset 4 
(3 meter level). 
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Eee Ee RIOD 1 

The highest mean wind speed of any period examined was 
observed during period 1. Observations were made at the 6 
Geamemereycl, Wave heights were between six and eight feet, 
with a mean wind speed Hanus coews Tits perred 1s, eharac- 
domeeed DY a complete departure from potential flow predictions. 

iitesphase-amplitude results for u, though basically 
Sinusoidal in subsets 1 and 2, Figures 17 and 18, show a 
Mieeer shitt back over the upwind side of the wave crest 
Gieiad Deen observed in period 16. The shift is sti1] 
Pueaomie in Subset 3 (Figure 19), and to a lesser extent in 
E@ieet 4 (Figure 20); however, the shape of the wave-form 
has become non-sinusoidal. The amplitude of the horizontal 
fluctuations are observed to be relatively small. 

A non-negligible role for turbulence is evident throughout 
Gmempertod for w,. the wave shape of the vertical velocity 
iPetecwatlOns varies in all subsets of the period and 
facooraal forms are not observed. The amplitude of the w 
component is athaceneey diminished from all other periods examined. 

The wave-related stress (uw) varies considerably over the 
period, ena changes of up to 90° are observed from subset 
Meeomsubpset 4, 

This distinct departure from potential flow, indicated by 
the non-sinusoidal shapes observed in the results of this period, 


could be due to non-linear processes. 
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Figure 17. Phase-Amplitude Results for Period lL zubset)).. 
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Figure 18. 
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Phase-Amplitude Results for 
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Period 1, Subset 2Z. 








SMALL LARGE TOTAL 


Pigure 19. Phase-Amplitude Results for Period 1, Subset 3. 


46 








SMALL LARGE TOTAL 


Figure 20. Phase-Amplitude Results for Period 1, Subset 4. 
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See EERIOD 6 

Meestirements in perirods6, were made at 8 meters and it is 
emeracterized by the low mean wind speed of 7.0 M/sec and 
the lowest observed wave heights at three to five feet. 

Wave heights had been diminishing during the previous six 
Memmes and wind speeds had been increasing slightly up to the 
Meme of the period, then diminishing. 

Bile=siountreucant teature Of this period, as indicated in 
Heeores Zi, 22, and 23 1s that of all periods examined it 
Seems to be very nearly as predicted by potential flow 
Gieory., ihis is especially discernable for the total waves 
phase-amplitude. However, there are some departures from 
potential flow in the phase relationships of u and w and in 
Ewe amplitudes of uw. 

Gnessivelht Shift torward of the negative extrema of the 
W fale tuations will also be noted in period 14 and will be 
Meeeriped in detail in the discussion of that period. 

Ime subset 1 (Figure 21) there is a shift forward in the 
ime cast 1ve extrema of the u component while the positive 
Ceenema 1S Consistent with results of period 12. The positive/ 
negative extrema INgetoOr the vretlarnder of the period are~ 
peso consistent with previous discussions in period 12 
mmeene o meter level, 

Potential flow theory predicts equal magnitudes (Figure 1) 
of uw in all quadrants. In these results the amplitudes 


ie sOveretine Wlase. of the wave and from subset 1 to subset 3. 
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In subset 1, Figure 21, quadrant IV, the amplitude ace 
are greater than the other quadrants as evidenced by the 
M@eecominance of the downward transfer of excess momentum. 
Peaeene end of the period the dominant quadrant is quadrant II, 
migbertive Of an Upward transfer of deficit momentum. In 
MemeEOtal wave, however, quadrant II dominates in subsets 1 
pimemez, Shifting to quadrant IV in subset 3. The total wave 
uw Closely resembles that predicted by potential flow theory. 
Nee aeterminations concerning FLIP's motion seem possible 
from this period; however, with the small observed wave 
Nerehts and a minimum of wind it is considered that the 


Pimors introduced were minimal. 
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Phase-Amplitude Results for Period 6, Subset l. 


Preaure 21. 
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LARGE 


SMALL 


Enase-Ampiitude Results for Period 6, Subset 2. 


meeoure 22. 
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Fasure 23. Phase-Amplitude Results for Period 6, Subset 3. 
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eee PERIOD 14 

ies perlod 1S unique in@that@neasurements weme made at 
the 3 meter level and thus offers a comparison with the lower 
level of period 12 and Thompson's (1972) results of the same 
data. 

meni ficant features include distinct deviations from 
Howential flow theory in the phase of w and the phase and 
amelatude of uw. Features of this period also deviate 
eemoracrably from those for the 3 meter level in period 12. 
However, in comparison to period 12, this period's mean wind 
Speed was only 7.6 M/sec and the observed wave heights were 
Gapeesix to eight feet. 

The phase relationships of u are comparable to potential 
imewerim Subsets 1, 2, and 3 (Figures 24, 25, and 26), 
emeeme for the consistent shift to the back node of the crest 
of the wave. However, in Figure 27, subset 4, it can be seen 
Mmeeemene Shitt Of U 1S to the front side of the crest. 

iiompson had noted the occurrence im this period of a 
minimum value of w before the inflection point of the wave. 
However, with the increase in resolution, this shift is 
Cons raered minor and Should afurther anerease in resolution 
be achieved it is doubtful that the shift would be discern- 
mec at ali, There 1s a considerable shift forward of the 
positive extreme of w which is a distinct departure from 
mecential flow predictions where the maximum w values should 


eccur at the inflection point of the wave at the end of 


Si 





aiia@rant IV. The shift observed in this period varies 
from 30° to 90° forward and accounts for the predominance 
femene downward transfer of excess momentum in quadrant IV. 

The discussion by Thompson (1972) of the uw phase- 
paiplertude results is consistent with these results from the 
same data and will not be repeated here. 

The differences between this period and the 3 meter 
Baeereor period 12 seem to be the result of the reduced wind 
freeads and wave heights. In the u component there is no 
mee rence between the phase relationships of the two periods. 
iiesanplitude of the u component fluctuations is greater in 
the period with the lower mean wind speeds and wave heights, 
fimeeiemay also explain the differences in the phase of w 
Memmeen the two periods, It does not seem likely that FLIP's 
Nemton would have led to these differences. The results in 
period 12 are considered to be consistent with measurements 
meomed £1xed platform by Davidson and Frank (1973) which 
fememalso in agreement with recent turbulence models. If 
Pond's (1967) suggestion concerning the tilt of the sensor 
Mmltenmorm 1S considered, the reduced tilt associated with 
meduced wave heights would certainly not lead to increases 
mimene OoScillations of FLIP and introduce an increased error. 
These differences in the results must then be attributed to 
Bmeilanging Significance for the role of turbulence at various 


combinations of wind speed and wave height. 
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Frpure 24. Phase-Amplitude Results for Period 14, Subset 1. 
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Figure 25. Phase-Amplitude Results for Period 14, Subset 2. 
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Phase-Amplitude Results for Period 14, Subset 5. 


Peoure 206. 
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Phase-Amplitude Results for Period 14, Subset 4. 


Feoure 27. 
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VI. SUMMARY AND CONCLUSIONS 


Analyses were performed on observational data obtained 
mes O©CeCan Waves from a relatively stable floating platform. 
The analyses on these data provided results on the phase and 
amplitude relationships of wave-related velocity fluctuations 
and wave heights. The results were interpreted with respect 
mompoetcential flow predictions for a wave propagating beneath 
fone ar tlow and the role of turbulence in wind-wave coupling. 
Interpretations were also made with respect to the possible 
inmerucnce FLIP may have had on the data and the results of 
the analysis of that data. 

Peesienificant feature. was the consistent phase relation- 
eimpeo: the velocity component u, with that of the wave. In 
elm@st all cases there was a shift of 15° to 30° of the nega- 
tive/positive extrema of the u component to the upwind side 
@emene wave crest/trough, Although this feature is very 
eeeelar to that predicted by potential flow theory, it is 
considered to be a significant deviation from the 180° pha e 
Giterence predicted. Also observed was the absence of the 
a@adrature relationship of u and w as predicted by potential 
miew, Deviations from potential flow predictions were 
interpreted to be indicative of interaction between the wave- 
induced motion and the shear flow. Thus a provision for 
mimbpulence interaction should be considered essential to 
Wind-wave coupling models. 


In comparison to potential flow predictions which have 


equal magnitudes of uw in all quadrants, the uw phase-amplitude 


Sie. 








fmooults revealed variations in magnitude over different 
phases of the wave. This suggests that nonlinear processes 
are responsible for momentum transfer in wind-wave coupling. 
Brom Comparisons with Similar results from a fixed 
sensor platform and an examination of the phase relationships 
Mmerne Dresent results, it was concluded that the floating 
mreerument platform had no appreciable effect on these 
results. However, results from similar investigations with 
heal time measurements of platform motion are needed to 


prrcemethen thas conclusion. 
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A eed ei BNO Ar 


ilps Motion will be described relative to a fixed 
Mmemenence coordinate system which has its origin at the 
equalibrium point of the buoy's metacenter (Figure 28). The 
moms pitch, roll, and yaw are the angular displacements 
Mectations) about the x, y, and z axes, respectively. Heave- 
imemon refers to displacements in the vertical (z) plane. 
immo these definitions the velocity of any point on the 


Mmievyecan be found relative to the fixed system, 


(1) Uses 0!) + ure a4 


mice Tc Up, Ese ties resultant Vector cr Ocean. v4 DOI. Onwere 
boys due to the linear velocity, Ur, of the buoy with respect 
fomctne fixed origin, and Ox T (Figure Z8) is the tangential 
momocity of the point about this origin, Measurements are 
mien with respect to the moving system, so the apparent 
welocity relative to the moving point on the buoy is 


(eZ) U' =U, -oxr 
R 


This is the velocity detected by the anemometer on board 


fie. if there is no wind. If a "real" wind oe ried 


1 appendix A is unpublished work by Dr. A. W. Green, Department 
Stewfeteorology and Oceanography, University of Michigan. It 
Weaeebeen revised slightly for inelusion in this thesis. 
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Figure 28. The fixed reference coordinate system has its 
origin at the metacenter of the spar buoy. For 
small pitch-roll displacements the metacenter 
Mommuced (Sec 1: Mand sa stages. p.mac., 
Canbridce ress. 19605) 
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Moving past the buoy is added, and only apparent velocities 
due to the heave (U,, = kw, ) plus rotation about the meta- 


Semeer are considered, thin (2) may be rewritten as 


(2a) Unseen, “Gera + ki = wo x 1 


where Wh is the oscillation along the vertical due to heave 
a ° * 

femecd by the anemometers. be is equal to the instantaneous 

mieal'’ wind vector. Now equation (2) is placed in component 


form . 


(a) u' = uae (w 52 - Wy) 


(i) Ci = = (w 2X - Ww 2) 


(c) w' = CS wie (wiy - WX) 


Site quantities Wy» W Ww, can be defined as the roll, pitch, 


>? 
and yaw rates; x, y, z are the instantaneous distances along 
jMmiemrespective axes for the equilibrium metacenter. The 

drift of FLIP due to the currents and wind will be ignored, 
eeenough it could be important in cases of very low wind 
fimecities. Next consider some of the statistical properties 
emene apparent wind field. 

The power spectral density is usually computed from the 
pea- correlation of a time series of the velocity fluctuations, 
which may be represented by 
p72 


(4) Ru! (t) = lim, A ff ourct) ut(t) at 
| -T/2 
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where Bee, Cop) ES stUtO-CommciatlOny HmeCluLOmTO mate apparent 
x-component of the measured wind field with the time lag, 


ijeeeine power spectral density is 


(5) G (f) = 2 f R 


Semerent sinusoidal motions will then give peaks in the 
power spectral density plots. Taking two examples, the 
effects of buoy motion on the power spectrum will be shown. 
ieee iic auto-correlation function of the horizontal | 
PUG tuat ON Component at zero lag (7 = 0) (the 


Paweimdicates time averace) . 





Pans iB <n 2 eee 
(6) Ra (tT=0) = utul = CO Et du (w,2-wzy) a (w52 W2Y) 


(1) (2) (3) 


ie across CoOtrelation. 2 uncti on. On thewnoricomea | 


ana verticals Fluctuations. .0O8 Reynolds Stress. 





coe’ = Uw + UW) - U (wWyy-w,S) - W (w52-wy) 
(1) (2) (3) (4) 


“Wy (wY-W5X) + (w2-wzy) , (WY -W5X) 
(5) (6) 
Mme first term on the r.h.s. of (6) is the mean square 


menctuiation, or the ''true zero-lag."' The second is the 


correlation of the wind fluctuation with the pitch and yaw 
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Memuons, which could have the effect of decreasing the 
meparent mean Square values. Finally the last term repre- 
Semmes the apparent net mean square of the correlated motion 
of pitch and yaw, which must be positive. It should be 
Memed that the yaw and pitch velocities are the total 
Mictamtaneous motions due to the rotation of the sensor 
peeme about the reference axes. Consequently, the terms 
buoy motion represents are the instantaneous apparent velocity 
components which by themselves have a virtually continuous 
Becetrum. According to Rudnick's (1967) measurements of 
ees MOtLon, the spectra tend to have peaks which roughly 
eemrespond to the wave frequency, pitch-roll, and heave 
resonances. Also the nonlinear response of the buoy can 
leads to higher order responses which would give peaks in power 
cameeveral Combination frequencies of the above. 

Mien rrrst term on the rih.s. of (7) ats: the "real" 
Peympodds stress. the second term could be troublesome si <e 
MPmemapsolute magnitude could be large, if there is a wind 
component which correlates Significantly with the vertical 
femeeity due to heave. The third term is the correlation of 
the wind fluctuations ance nenVetiler EevVCloOGitdes sdile. CO“piseen 
and roll. The fourth is the correlation of the vertical wind 
matetuation and the pitch and yaw motions. The fifth term 
memchne correlation of yaw-pitch with roll-pitch motions of 
miempuoy, thrs term is probably negative since the pitch- 


perehwcorre lation, (w,2) ° (-wx) is less than zero. BOMEX 
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Gera MOtilon Spectra have relatively sharp peaks at lower 
fPeeguencieés With lower values for u'w' than other observers. 
Empeobable reason is the presence of the fourth and fifth 
terms of (7) which may have negative values that tend to 
decrease the mean Reynolds stress. 

moparent FLIP motions, which are more or less sinusoidal, 
iyeope conveniently represented as the super-positions of 


mitemeomponents of the apparent velocity in the following way, 


E{ (iw 


7 ee x [iX Cos (w 


lente 1, nt? 


™ 
e 


t 
nt . - ¥ So ene * ¢ n? 


io 


Kon nnSINGy nt * ony)? 


winere SUMNMation 18 over all modes of FLIP. 


= are random phase variations; a a) OO ae ea ea 


t 
n> Ons Onn BE | 
the maximum displacements from equilibrium position associated 
With a given mode. 
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